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Abstract: The scattering theory of Lax and Phillips, originally developed to describe res- 
onances associated with classical wave equations, has been recently extended to apply as 
well to the case of the Schrodinger equation in the case that the wave operators for the 
corresponding Lax-Phillips theory exist. It is known that the bound state levels of an atom 
become resonances (spectral enhancements) in the continuum in the presence of an elec- 
tric field (on all space) in the quantum mechanical Hilbert space. Such resonances appear 
as states in the extended Lax-Phillips Hilbert space. We show that for a simple version 
of the Stark effect, these states can be explicitly computed, and exhibit the (necessarily) 
semigroup property of decay in time. The widths and location of the resonances are those 
given by the poles of the resolvent of the standard quantum mechanical form. 
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1. Introduction 

The description commonly used for an unstable system is that of Wigner and 
Weisskopf 1 . They assumed that the unstable system is represented by a vector, say <fi, 
in the Hilbert space of states, and that the Hamiltonian evolution of this state permits the 
development of other components which represent the "decayed" system. One may assume 
a Hamiltonian Hq for which <p is an eigenstate, and the full Hamiltonian is constructed by 
adding a pertubation V which induces transitions from this state. The basic assumption 
of their physical model is that these transitions carry the state of the system from that of 
some specific system to a state in the continuous spectrum which contains the decay prod- 
ucts. For example, one may think of an atom in an excited state as the unstable system; 
in the absence of electromagnetic interaction, this state would be a stable bound state. 
A term added to the Hamiltonian corresponding to electromagnetic interaction induces 
a transition from this state, and the square of the corresponding transition amplitudes, 
which may be generally computed perturbatively, give the probability for decay. This very 
standard technique is rigorously correct for reversible quantum transitions, according to 
the laws of quantum theory. When applied to the decay of an unstable system, however, 
for which the evolution is irreversible, it consitutes an approximation which may be in- 
adequate. In the following we discuss some of the difficulties of the application of the 
Wigner- Weisskopf method to the treatment of irreversible phenomena, such as particle 
decay, and in Section 3 we describe the general structure of the Lax-Phillips theory, which 
yields an exact semigroup evolution and appears to be a more appropriate description of 
such phenomena. In Section 2, we apply the Wigner- Weisskopf theory to a simple model 
for the Stark effect, for comparison, and work out the Lax-Phillips theory for this model 
in Section 4. Some of the calculations in Section 2 coincide with those needed in Section 4; 
the approximate pole approximation decay law of the Wigner- Weisskopf theory contains 
exactly the same pole as the singularity of the Lax-Phillips S matrix, so that the decay 
laws are (approximately) the same, but, as we shall see, there are important differences in 
that the Lax-Phillips treatment leads to an exact semigroup and contains information not 
accessible in the Wigner- Weisskopf method. 

One easily sees that the Wigner- Weisskopf procedure cannot yield a semigroup prop- 
erty of evolution. The "survival amplitude" for such a state is given by 

A{t) = {^e- im ^ (1.1) 

and the survival probability is |A(t)| 2 . One observes, in general, that decaying systems 
are characterized by exponential decay in time. However, the time derivative of |A(t)| 2 at 
t = is 

j t \A(t)\%=o = -ifaHfi + HHM) 
= 0, 

where the last equality follows from the self-adjointness of H. Hence the probability curve 
must arrive at the origin with zero slope, and cannot be exponential. Generally, this 
deviation from exponential occurs only for very short times 2 ' 3 , e.g., for a weak decay with 
lifetime of the order of 10~ 10 — 10~ 8 seconds, this short time may be of order 10~ 22 seconds. 
The dominating exponential behavior for intermediate times can well account for the data 
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in many such cases of single channel decay. For multichannel problems (many types of 
final states), for example, the two channel neutral K-meson decay, the lack of semigroup 
property may lead to disagreement with experiment. 

To see this (and to see as well the origin of the exponential part of the decay), let us 
examine briefly the two channel case. In this case, the probability of decay is given by 

P = J2 I ^|<A,|e-^V)| 2 , (1-3) 

3 

where {Xj} corresponds to the continuous spectrum of the "final" states, and <fi is an 

arbitrary state in the two dimensional subspace of initial states (e.g., K° 7 K°). This can 
be written alternatively, since the evolution is unitary, as 

P= 1-£K^' e- iH V)| 2 , 

3 

where {(f) j} spans the subspace. Since one can always express as a linear superposition 
of the {4>j}, it is sufficient to study the quantities 




where W(z) is a 2 x 2 matrix defined by the second equality, C is a contour for the inverse 
Laplace transform, and the quantity (z — H)~ x is the resolvent operator associated with 
the Hamiltonian H (defined for z in the upper half plane). The contour C runs above the 
real line from +oo to — oo. Almost every 2x2 matrix has a decomposition 4 of spectral 
type 

W(z) = z 1 (z)Q 1 (z) + z 2 (z)Q 2 (z), (1.5) 

where the Qi(z) are constructed from the right and left eigenvectors of W(z), and satisfy 
Qi\ — Qh QiQj = 0,i ^ j. The resolvent kernel is regular in the upper half plane, but 
its analytic continuation through the cut induced by the continuous spectrum of H may 
have poles (or, in principle, other singularities). For small perturbations, one expects that 
the bound states in Hq for a model for which H = Hq + V would move slightly into this 
region, which we shall call the second Riemann sheet 4 . Poles would appear as zeros in 
either z — z\{z) or in 2 — 2^2(2;), at points z a or Zb- Let us assume that the eigenvalues 
are ordered so that one pole appears at z a such that z a — zi(z a ) = and the second at Zb 
such that Zb — Zi(zb) = 0. In distorting the contour below the real axis across the cut into 
the second sheet, these pole contributions may dominate the intermediate time behavior 
of the amplitude (the very short time behavior is controlled by the Taylor expansion at 
small t) s . For semibounded spectrum, in (0, 00), the long time behavior is controlled by 
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the contribution at the branch point, and goes as t n , where n is the dimensionality of 
space. The resulting pole approximation is then 



A l3 (t)=e 



—iz a t 



Qi(z a ) + e 



-iz b t 



Q2(z b ), 



(1.6) 



where we have also approximated the factor in the residue arising from the derivative of 
the denominator in the neighborhood of the singularity by unity. 

For the single channel case, one sees the emergence of the simple exponential behavior 
for intermediate times. However, for the two channel (or more) case, we see that the matrix 
Aij(t) does not satisfy the semigroup property 



because, in general (in the absence of CP conservation), Qi(z a )Q2(zb) ^ (in case W(z) 
is independent of z, these matrices would be, of course, orthogonal). 

One may think of this semigroup property as a requirement for the evolution of an 
irreversible system 5 ; this follows from arguing that the evolution on a Hilbert space must 
be composable (in the sense that a product of evolutions must be equivalent to a single 
evolution over the total interval) for a sequence of i's that are positive, since the system 
carries no dependence on its past, but the inverse does not exist. Indeed, experiments 
carried out at Fermilab 6 show that the phenomenological model of Lee, Oehme, Yang 
and Wu 7 is an extremely accurate representation. This model assumes that the evolution 
in the two dimensional subspace is given by exp(— iH2x2t), where #2x2 is a 2 x 2 non- 
Hermitian matrix. In this case, the matrix H 2X 2 , independent of z, replaces W(z) in 
the the spectral decomposition (1.5); the pole residues are orthogonal, and one obtains 
the exact semigroup law. Using a Lee-Friedrichs model 4 to estimate the errors, one finds 
that the non-orthogonality of Q(z a ) and Q(zf,) induces errors that are larger than the 
experimental errors in regeneration experiments and in the interference in the exit beam 8 . 
We infer from these results that for irreversible transitions (of the type occurring in particle 
decay), one should have the semigroup property. 

The theory of scattering developed in 1967 by Lax and Phillips 9 for application to 
classical wave equations has the property that resonance decay satisfies exactly the semi- 
group property. It is also true that, in that theory, there is a well-defined state in the 
Hilbert space corresponding to the resonance. 

This theory has been developed and methods found to apply it to resonances in the 
quantum theory 10-14 . In the Wigner-Weisskopf method, the resonance is represented only 
as a pole in the complex plane. Attempts have been made to extend the Hilbert space to the 
so-called rigged Hilbert space 15 (Gel'fand triple constuction) , in which exact exponential 
decay occurs, but the "eigenstate" of the resonance, in that case, lies in a Banach space 
that is not generally a Hilbert space, so that scalar products and expectation values are 
not defined. The definition of the wave function of the resonance carrying information on 
its size, shape and momentum distribution could be important in many applications, for 
example, for the resonance in a quantum dot or a small metallic particle 16 . 

It is most straightforward to apply these methods to problems in which the correspond- 
ing Hamiltonian of the standard theory has spectrum — oo to +oo. One such example is 



A{t) A{t') = A(t + t') 



(1.7) 
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the relativistic quantum theory 7 , where the unperturbed Hamiltonian has a form pro- 
portional to E 2 — p 2 ; the Lax-Phillips theory for such a problem has been worked out 18 . 
For the non-relativistic case, the spectrum of the Hamiltonian is generally bounded from 
below, and further studies are being carried out in order to be able to apply the theory 
in these cases as well 14 . However, the non-relativistic Stark effect is a case in which the 
non-relativistic problem has a Hamiltonian unbounded from below. We shall study this 
problem in detail in this paper, and show that one finds semigroup behavior for the de- 
cay of the Stark induced resonances; the wave functions for the resonant states are also 
explicity worked out. In the following we study a simple model, of the type of Friedrichs 
and Rejto 19 , for the Stark effect, in the framework of Wigner and Weisskopf, and then 
imbed the results into the more general framework of Lax and Phillips. We shall review 
the content of the Lax Phillips theory in this later section. 



5 



2. Wigner-Weisskopf Analysis of a Stark Model 

The potential for the Stark effect problem, of the form —Ex, is unbounded (on the 
full space; we shall work in one space dimension). For a model of the form 19 

H = -Ex + AP , (2.1) 

where A is real and Pq is a rank one projection operator, it will be convenient to con- 
sider —Ex = Hq as the unperturbed Hamiltonian and the second term, APo = V as the 
perturbation (the resolvent is, of course unaffected by this choice, but the form of the 
perturbation theory is very different) . We study in this section the Wigner-Weisskopf de- 
scription of the resonance, and in the next section, imbed this analysis in the Lax-Phillips 
Hilbert space. 

Let us choose for P the form 

<x\P \x'>= {-y e - {x2+x ^ (2.2) 

TV 

The resolvent satisfies the identity (second resolvent equation) 

G = G + G VG, (2.3) 

where, as above, G = (z — H)~ x , and Go = (z — i^o) -1 , defined for z in the upper half 
plane, where Hq = —Ex. The x, x' matrix element of G is therefore 

1 1 r°° 

(x\G\x') = — 5(x - x) + —A / (x\P \x")(x"\G\x')dx". (2.4) 

z + Ex z + Ex J _ oc 

Let us define 

f(z, x') = J dx"e- x " 2 (x"\G\x')dx". (2.5) 
It then follows from (2.4) that 

2 2 I 

f(z,x') = [ dx^-— + \ I dx^—J-e- (x2+x " 2 \x"\G\x')dx". 

J J Z + EX J Z + Ex\l 7T 



One can write this as 



or 



where 



e x 



/(*. x ') = 7r VTW^ (2 - 6) 



-2t 2 

e ITT _ 2z£ „ r . /- Z 



F(z) = j dx^^ = -e-^erfc[iV2-}. (2.7) 



Returning to Eq. (2.4), we see that 



1 e , /N , /2 e" x e 
-5(x - x ) + A 



-x 2 -x' 2 -i (2.8) 



^ + -Ex V 7T (^ + £ , x)(z + Ex') 1 _ ^ f2.F{z) 



We now wish to approximate the time behavior of the survival amplitude. As in 
Eq.(1.4), the time dependence of the survival amplitude (for one channel) is given by 

A(t) = ^-j^\G(z)\^e-- t . (2.9) 

The contour C corresponds to a line running in the complex energy plane from right to 
left slightly above the real axis. The matrix element (ip\G(z)\ip) is analytic in the upper 
half-plane. We can shift this line continuously and differentiably through the real axis into 
the lower half plane, provided that the contribution of the vertical pieces at ±oo vanish. 
It is clear from (2.7) and (2.8) that this is true for the part of the vertical integrations that 
lie in the upper half plane. To write the integrand along the new curve below the axis, we 
must analytically continue F(z). To do this, consider (for £ real) 

/oo 11 
dxe- 2x2 {- — — - 



= -2ni J dxe~ 2x2 5(£ + Ex) 



(2.10) 



2ni a gf 
-e 



E 

This function has an analytic extension in the finite lower half plane, given by 

F\z) = F(z)- 2 -^e-^. (2.11) 

We find, numerically, for a reasonable choice of parameters and a simple assumption 
for y?(x), that the analytic continuation of the function 



((p\G(z)\(p) = I dxdx'cp*{x)( — l—S(x-x') 
J V z -\- Ex 



x 2 „-x' 2 i v (2.12) 



2 e~ x e~ x 1 
+ A V * (z + Ex){z + Ex') " _ x fI F(z) ^ 



defined by (2.8) and (2.11), into the lower half plane, has a pole, inducing an exponential 
decay term to the amplitude. The function (<p\G(z)\<p) has a very simple form if we assume 
that (p(x) has the Gaussian form 

V?(x) = yV*\ (2.13) 
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The first term of (2.12) contains e~ 2x ; its integral with the denominator z + Ex is, 
according to (2.7), the function F(z). The second term factorizes into two integrals of the 
same form. It then follows, with this assumption on <p, that 



The analytic continuation of this function into the lower half plane is achieved by the 
continuation of F(z) to F (z); this function has no poles in the finite lower half plane, and 
hence the pole can only come from the condition 



For the value X/E = 11, MAPLE provides us with a unique solution for the position 
of the pole, zq = —4.446 — .31896 x 10 _15 z, which has, as expected, a very small imaginary 
part for this reasonably physical choice of parameters. Differentiating (2.15) implicitly 
with respect to E, we find that the real part of the pole moves to more negative values as 
E increases. Since the unpertubed system has mean position of the particle at zero, this 
shift corresponds to an increase in field induced polarization with increasing value of the 



There remains, however, a contribution to the survival amplitude from integration on 
a line running from +oo to — oo on the real part of z = £ + i£, where £ can be very large 
and negative. The contribution of this so-called background integral is strongly suppressed 
by the exponent exp(— izt) for t large and positive. For small t, however, this suppression 
is not strong unless £ — > — oo. However, in this limit, the integrand is not well-defined, 
since for any large and negative £, the contribution of the discontinuity in F(z) strongly 
suppresses the integrand for £ small compared to (. This suppression is not maintained, 
however, for the contributions from £ in the neighborhood or greater than (. Hence the 
convergence of the contribution on the background is not uniform. We see from the general 
argument (1.2) that the pole contribution cannot represent the result precisely for small t, 
and therefore the Wigner-Weisskopf treatment, even in this case of unbounded spectrum, 
cannot result in a pure exponential (semigroup) behavior for the reduced evolution. It is 
exactly in this respect that the Lax-Phillips treatment provides a result which is closer to 
the physics of irreversible decay. 

3. Lax-Phillips theory 

The Stark effect admits a simple transition to the Lax-Phillips framework, since the 
Hamiltonian is unbounded. In the following, we briefly describe the theory of Lax and 
Phillips 9 and its extension to the quantum case 13 , and then apply the formalism to the 
model we are using. 




(2.14) 




(2.15) 



field. 
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The scattering theory of Lax and Phillips assumes the existence of a Hilbert space H 
of physical states in which there are two distinguished orthogonal subspaces T> + and £>_ 
with the properties 

U(t) V+C V+ t > 
U{t)V- C V- t < 

f]U(r)V ± = {0} (3 .i) 

T 

\JU(r)V± = H, 

T 

i.e., the subspaces T>± are stable under the action of the full unitary dynamical evolution 
U(r), a function of the physical laboratory time, for positive and negative times r respec- 
tively; over all r, the evolution operator generates a dense set in H from either V + or £>_. 
We shall call V + the outgoing subspace and £>_ the incoming subspace with respect to the 
group U (r). 

A theorem of Sinai 20 then assures that Ti can be represented as a family of Hilbert 
spaces obtained by foliating Ti along the real line, which we shall call {s}, in the form of 
a direct integral 

H= I H 8 , (3.2) 

J® 

where the set of auxiliary Hilbert spaces Ti s are all isomorphic. Representing these spaces 
in terms of square-integrable functions, we define the norm in the direct integral space (we 
use Lesbesgue measure) as 

/oo 
MMh, (3-3) 
-oo 

where / e H represents a vector in Ti in terms of the L 2 function space L 2 (— oo, oo, H), and 
f s G H, the L 2 function space representing H s for any s. The Sinai theorem furthermore 
asserts that there are representations for which the action of the full evolution group U(t) 
on L 2 (— oo, oo, H) is translation by r units. Given D± (the 1? spaces representing T>±), 
there is such a representation, called the incoming translation representation 1 , for which 
functions in D- have support in L 2 (— oo, 0, H), and another called the outgoing translation 
representation, for which functions in D + have support in L 2 (0, oo, H). 

Lax and Phillips 9 show that there are unitary operators W±, called wave operators, 
which map elements in H, respectively, to these representations. They define an S'-matrix, 

S = W+WI 1 (3.4) 

which connects these representations; it is unitary, commutes with translations, and maps 
L 2 (— oo,0) into itself*. The singularities of this S'-matrix, in what we shall define as the 

* Note that the S matrix of the Lax-Phillips theory maps the incoming to the outgoing 
representations unlike the S matrix of the usual scattering theory which maps incoming 
to outgoing states 
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spectral representation, correspond to the spectrum of the generator of the exact semigroup 
characterizing the evolution of the unstable system. 

With the assumptions stated above on the properties of the subspaces V + and £>_, 
Lax and Phillips 9 prove that the family of operators 

Z(t) = P+U(r)P. (r > 0), (3.5) 

where P± are projections into the orthogonal complements of T>±, respectively, is a con- 
tractive, continuous, semigroup satisfying 

Z(t 1 )Z(t 2 ) = Z(t 1 + t 2 ) (3.6) 

This operator annihilates vectors in T>± and carries the space 

ic = nev + ev_ (3.7) 

into itself, with norm tending to zero for every element in /C. One may demonstrate the 
semigroup property by noting that the operator acting on a vector in D_ vanishes by 
definition; on a vector in D + , it vanishes due to the invariance of D+. The scalar product 
on the left with a vector g in D + vanishes by definiton; if g is in D_, the factor P + may 
be suppressed, since such a g is in a subspace othogonal to D + . Hence the operator (3.5) 
can only connect /C to /C. In the product 

P + U(n)P-P+Ufo)P- (3.8) 

one may suppress the factor P_ following P + between the two unitaries since vectors in /C 
are orthogonal to _D_ . Then, consider 

P + U(n)P+U(T2) = P+U(n)[i - (l - P+)]U(r 2 ); 

the second term in square brackets makes no contribution since it is a projection into D+, 
and U(ti) leaves this subspace invariant. The composition of U{j\) and U{j2) then proves 
(3.6). 

The outgoing subspace D + is defined, in the outgoing representation, in terms of sup- 
port properties (this is also true for the incoming subspace in the incoming representation). 
One can then easily understand that the fundamental difference between Lax-Phillips the- 
ory and the standard quantum theory lies in this property. The subspace defining the 
unstable system in the standard theory is usually defined as the eigenstate of an unper- 
turbed Hamiltonian, and is not associated with an interval on a line. The subspaces of the 
Lax-Phillips theory are associated with intervals (i.e., the positive and negative half-lines 
in the outgoing and incoming representations). It was shown in ref. 11, for a Lax-Phillips 
theory constructed on a model of free evolution, that the generator of the motion re- 
stricted to the subspace /C is symmetric but not self-adjoint. To briefly review this point 
in a slightly different way, consider the infinitesimal generator of the semigroup. Accord- 
ing to the Stone theorem, the unitary evolution U{r) can be represented as e~ tKl ~; since 
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P + U(t)P- takes /C into itself, as shown by Lax and Phillips, P + U(t)P- = P)cU(r)P)c, 
and therefore the infinitesimal generator (the derivative at r = 0) satisfies 



P+KP_ = P K KP K . 

Consider vector-valued functions f(s),g(s) in /C, in the translation representation, so that 
(we may suppress _P_ acting on / or g) 

(P + K)f(s) = -i6(-s)^-f(s). (3.9) 
We then take the scalar product with g to obtain 



I g*(s)(P + Kf)(s)ds = -i [ g*(s)9(-s)- 

J — oo J 



ds 



ds 



= -ig*(0)f(0) + J {-i6(-s)^yf(s)ds 
= -i</(0)7(0) + J(P+Kg)*(s)f(s)ds, 



so that the operator P + KP- = PjqKPjq is not self-adjoint. It is through this mecha- 
nism, that the subspaces are defined by (semibounded) support properties on a variable 
acted upon by an operator represented as a derivative in the translation representation 
(somewhat analogous to the non-self-adjointness of the operator id/dr in three dimen- 
sional spherical coordinates), that the Lax-Phillips theory provides a description that has 
the semigroup property for the evolution of an unstable system. 

It follows immediately from the property of a contractive semigroup that the generator 
has, in fact, a family of complex eigenvalues in the lower half-plane; the eigenf unctions are 9 

where n is some vector in the auxiliary space. 

In order to construct an imbedding of Lax-Phillips theory into the quantum theory 13 
we assume that there exist wave operators 0± which intertwine this dynamical operator 
with an unperturbed dynamical operator Kq. We shall assume that Kq has absolutely 
continuous spectrum in (— oo, oo). 

We begin the development of the quantum Lax-Phillips theory with the construction 
of the incoming and outgoing translation representations. In this way, we shall construct 
explicitly the foliations required by the Lax-Phillips theory described above 13 . Let us 
define the free spectral representation 9 in terms of the Fourier transform 

f (a(3\g) = J s f (s(3\g)ds, (3.11) 

where f{a/3\g) satisfies 

f{c'p\Ko\9) = o'f{c-p\g)- (3-12) 
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Here \g) is an element of H and (3 corresponds to the variables (measure space) of the 
auxiliary space associated to each value of a, which, with a, comprise a complete spectral 
set. The functions may be thought of as a set of functions of the variables (3 indexed on the 
variable a in a continuous sequence of auxiliary Hilbert spaces isomorphic to H . Clearly, 
Kq acts as the generator of translations in the representation f(s/3\g). 

We now proceed to define the incoming and outgoing subspaces V±. 

Let us consider the sets of functions with support in L 2 (0, oo) and in L 2 (— oo, 0),and 
call these subspaces . The Fourier transform back to the free spectral representation 
provides the two sets of Hardy class functions 

f(vP\9o) = J ^ f(*0\9o)da G H ± , (3.13) 

for g± G Dq. _ 

We may now define the subspaces T>± in the Hilbert space of states H. To do this we 
first map these Hardy class functions in H to H, i.e., we define the subspaces Vq by 



/ 



Y,\<T0) f fWlgftdaeVf. (3.14) 



We shall assume that there are wave operators which intertwine K with the full 
evolution K, i.e., that the limits 

lim e iKr e- iKoT = Q± (3.15) 

T^±0O 

exist on a dense set in 7i. The conditions for the existence of such wave operators are 
satisifed if the difference between K and Kq is a "small" operator. We shall explicitly 
construct these wave operators in the case of the Stark effect that we study here. 

The construction of V± is then completed with the help of the wave operators. We 
define these subspaces by 

v + = n+v+ 



(3.16) 
V-=Q-Vq. 

We remark that these subspaces are not produced by the same unitary map. This procedure 
is necessary to realize the Lax-Phillips structure non-trivially; if a single unitary map were 
used, then there would exist a transformation into the space of functions on L 2 (— oo, oo, H) 
which has the property that all functions with support on the positive half-line represent 
elements of T> + , and all functions with support on the negative half-line represent elements 
of T>- in the same representation; the resulting Lax-Phillips S'-matrix would then be trivial. 
The requirement that V + and T>- be orthogonal is not an immediate consequence of our 
construction. Since the functions f{sj3\g^) have support in, respectively, the positive and 
negative half lines, and the orthogonality of T>± is determined by the integral of the product 
of these functions with an operator valued kernel S(s — s') (to be defined below), one 
sees that suitable analyticity properties of the transformed kernel S(a) assure that these 
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subspaces will be orthogonal. This analyticity property (upper half plane analyticity) is 
true in the Stark model that we treat. 

The wave operators defined by (3.15) intertwine K and Kq, i.e., 



KVL± = tt±K ; (3.17) 

we may therefore construct the outgoing (incoming) spectral representations from the free 
spectral representation. Since 

KV ± \a(3) f = n±K \a(3) f 
= an±\af3) f , 

we may identify 

\ap)out = n±\ap) f . (3.19) 

in 

The Lax-Phillips S-matrix is defined as the operator, on H, which carries the incoming to 
outgoing translation representations of the evolution operator K. Suppose g is an element 
of H; its incoming spectral representation, according to (3.19), is 

in((rp\g) = fiaplfiz'g). (3.20) 

Let us now act on this function with the Lax-Phillips S'-matrix in the free spectral repre- 
sentation, and require the result to be the outgoing representer of g: 

out {ap\g) = f {ap\n^g) 

= I 'da'J2f(°P\S\°'P')f fiv'P'^-Jg) (3 - 21) 

J f3' 

where S is the Lax-Phillips .S-operator (defined on H). Transforming the kernel to the free 
translation representation with the help of (3.11), i.e., 

f (s(3\S\s'(3') f = [ dada' e ias e- i,T ' s ' f {o-p\S\o-'P') f , (3.22) 

(27TJ J 

we see that the relation (3.21) becomes, after using Fourier transform in a similar way 
to transform the in and out spectral representations to the corresponding in and out 
translation representations, 



out(sP\g) = fism+'g) = I ds'YsfWWOff^'P'^g) 

0' 

(3.23) 

= / ds'J^fWlsWOfinWP'lg)- 

J at 



0' 

Hence the Lax-Phillips .S-matrix is given by 

S = { f (s[3\S\s'[3') f }, (3.24) 
13 



in free translation representation. It follows from the intertwining property (3.17) that 

f {af3\S\a'f3') f = 5(a - a')S^' (a), (3.25) 

This result can be expressed in terms of operators on H. Let 

wZ 1 = (3.26) 

be a map from 7i to H in the incoming translation representation, and, similarly, 

w^ = { f {sm- + 1 } (3-27) 

a map from 7i to H in the outgoing translation representation. It then follows from (3.23) 
that 

S = w~ 1 w_, (3.28) 

as a kernel on the free translation representation. This kernel is understood to operate on 
the representer of a vector g in the incoming representation and map it to the representer 
in the outgoing representation. 

4. Application to the Stark model 

Since the spectrum of the Stark model that we are using has spectrum (— oo, oo), we 
may take for the generator of motion 

K = Hstark] K =H Stark; V = Vstark, (4.1) 

where H stark , H stark and V = Vstark are the operators defined in Section 2. Since H 
is proportional to x, we may make use of the canonical commutation relations of the 
quantum theory to identify the momentum p as proportional to the (foliation) variable of 
the unperturbed (free) translation representation, i.e., one may take s = p/E, implying 
that [s, K ] = i (we have taken % = 1). We then have 

(e- iK ° T f)(s) = f(s-T), (4.2) 

or, in differential form, 

(K f)(s) = -i±f(s). (4.3) 

The auxiliary Hilbert spaces of the corresponding Lax-Phillips theory are one-dimensional. 

The spectrum of K Q , given by {—Ex}, with — oo < x < oo, can then be identified 
with a of the Lax-Phillips (unperturbed) energy representation. We shall follow this formal 
identification to develop the Lax-Phillips theory of resonances, and return to the original 
interpretation of x and p to obtain physical information about the resonant state. 

The wave operators are defined as 

Q± = lim e iKT e~ iKoT . (4.4) 

T^±00 
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We shall calculate the matrix elements of the wave operator in the unperturbed energy 
representation. It will be convenient, moreover, to use directly the measure on the spectrum 
of x; we therefore use kets of the form \x) = ^/p2j]<7). 

Following the standard procedure for taking these limits , one finds for the represen- 
tation that 

(x\n±\x') =S(x-x') - lim (x\ V\x'). (4.5) 

x 1 1 1 V ; e^o+ y l H + Ex'±ie 1 1 K 1 

Since this formula is bilinear in the kets \x) and \x'), one could use equally well the kets 
| a) and \a'). 

The operator multiplying V in (4.5) is —G(z), for z = —Ex'^fie; the matrix elements 
of this operator were evaluated in (2.8). Carrying out the integral for the product G(z)V 
with the help of the definition (2.7), one finds that 

(x\G(z)V\xf) = xJ = (4.6) 

^(z + Ex)(l-X^F(z)) 

The wave operators are then given by 

[2 e -(x 2 +x' 2 ) 

(x\fl±\x') = 8(x - x') + AW = (4.7) 

V 71 (E(x - x') =f - X^F(-Ex' t ie)) 

Using a partial fraction decomposition for the product of the GV terms, one easily verifies 
that the operators 0± are unitary. 

We now turn to the construction of the incoming and outgoing translation represen- 
tations. To do this, we follow the method of reference 13. We define the free outgoing 
translation as the set of functions with support in s (i.e., p/E in the Stark model) on the 
positive real axis. By (3.11), the functions 

/•OO 

f+(x)= / e*»/J(p)dp=(x|/J> (4.8) 
Jo 

are in the free outgoing translation representation, and are analytic in the upper half x- 
plane (lower half cr-plane). Since the wave operators intertwine Kq and K, functions of 
the (full) outgoing representation are then given by 

/oo 
(x'\n + \x)dxf°(x)dx = f° + u \x') E D+. (4.9) 
-oo 

Given a function of the type f+(x), one can calculate the resulting function f^! lt {x') 
explicitly by noting that the boundary value of F(z) from below the real axis is 

r°° e~ 2x " 2 in n r- 
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For the construction of the incoming translation representation, one uses fi_ and a 
corresponding set of functions /_ (x) with support on the negative half line. The kernel of 
integration then contains F(—Ex' — is); this may be obtained from (2.10), i.e., 

F(-Ex') above = F(-Ex') below - ^ e - 2x '\ (4.11) 

We now turn to the calculation of the S'-matrix. We see from (3.21) that 

(x\S\x') = (x\n+~ 1 tt-\x'}. (4.12) 

We now use the definition (3.15) for the wave operators; following the standard method 22 , 
we find that 

2tt? 

(x\S\x') = S(x - x')(l + — lim (x\T(-Ex' + ie)\x'), (4.13) 

E e^0+ 

where T(z) = V(l + G(z)V). We now compute 

lim (x\T(-Ex' + is) \x') = lim / \J - e -* 2 +*" 2 \5{x" - x') 

e^0+ s^0+J_ oo V 7T 



\Jl e -<- x +x ) 

-]dx" 



(E(x' -x")-ie)(l- X^F(-Ex' + is) ' ( 414 ) 

\ x fl e -(x 2 +*' 2 ) 

lim 

e->0_i 



+ 1- \J±F(-Ex' + ie) 

V 7T v ' 



It then follows that 



(x\S\x') = 8(x - x') \l + =^r I™ 1 , (4-15) 

E e ^ 0+ l-\JlF(-Ex + ie) 



so that we may write 

(x\S\x') = 8{x-x')S{x). (4.16) 

If we write the semigroup evolution, restricted to the subspace /C, as 

Z( T ) = e~ lB \ (4.17) 

it follows from the contractive semigroup property that the operator B has an eigenvector 
in the outgoing representation satisfying 

Bfout = H f out, (4-18) 
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with [i in the lower half plane, for which the eigenfunctions are of the form (with support 
in (0, oo)) 

fout(s) = e~^ s n, (4.19) 

where n is a vector in the auxiliary space (in our Stark model, one dimensional). The 
eigenfunctions in the outgoing x representation (corresponding to the "energy" variable 
a), are then of the form 

in 

fout(x) = • (4.20) 

X — fJL 

The S'-matrix, connecting the incoming to outgoing representations, therefore has the form 

S(x) — , (4.21) 

X — Zo 

where zq is the position of the pole of the diagonal .S-matrix in the lower half plane, which 
we identify with the semigroup exponent fj,, and r is the residue. From (4.15), one sees 
that the pole of the S'-matrix corresponds to a zero of the denominator 

1 - X^F(-Ex + ie) 
continued to the lower half plane, i.e., we must find the zero of 

l-\^-F l (z). (4.22) 

This is precisely the pole (zero of (2.15)) which occurs in the pole approximation of the 
Wigner-Weisskopf theory discussed in Section 2. The residue of the pole in (4.21) is then 
given by 

r~ " 2z o 



5. Conclusions and Discussion 

We have studied a model for resonances which provides the possibility of analysis 
using the Wigner-Weisskopf method, the standard technique for such analyses, giving 
exponential behavior in the pole approximation, but also accessible to direct analysis in 
the Lax-Phillips framework. In the Wigner-Weisskopf analysis, the resonance is described 
by the position of a pole in the complex energy plane, and does not have a state in 
the Hilbert space associated with it. Although there has been considerable study of the 
application of the method of rigged Hilbert spaces 23 , or Gel'fand triples (D C H C H, 
where H is defined as the set of bounded linear functionals on the subspace D), for the 
description of resonances which satisfies the property of exact exponential decay (without 
the "background" corrections to the Wigner-Weisskopf pole approximation), the elements 
of the Gel'fand triple are not, in general, vectors of a Hilbert space (they belong to a 
Banach space), and have no scalar products. Hence, it is not possible to compute the 
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expectation value of an observable (unless its extension to H maps it into the smaller 
space D on which the linear functionals of H are defined), or to study physical properties 
such as localization of the state. 

The Lax-Phillips formulation describes the resonant state as an element of a Hilbert 
space, and answers to such questions then become accessible. As a simple example, we 
have computed the resonant state for the Stark model that we have used in (4.20). The 
variable x used here corresponds to the "energy" in the Lax-Phillips formal structure, 
but retains its physical meaning in the result as position. Note that this is also true in 
our formulation of the Wigner-Weisskopf model, taking for the unperturbed Hamiltonian 
the term —Ex in the Hamiltonian (2.1) since this term is large compared to the term 
producing the imbedded bound state. The position variable occurs in the Hamiltonian 
producing a continuous energy spectrum {Ex}. The interpretation of the poles of the S- 
matrix, or resolvent therefore remains, as in the usual formulation of resonance problems, 
as occurring in the complex energy plane, but the variable x retains its physical meaning 
as coordinate as well. This method leads us to an approximate exponential decay law in 
time in the pole approximation of Wigner-Weisskopf theory, and to an exact exponential 
decay law in the Lax-Phillips treatment, with precisely the same exponent. 

In the framework of general Lax-Phillips theory, the resonant state carries the pole 
as in (4.20) where the x that appears there would be replaced by another symbol, say, cr, 
associated with unperturbed energy, and the distribution over space of the wave function 
would reside in the vector n of the auxiliary Hilbert space. In our case, this vector is just 
a number (one-dimensional), and the space distribution is provided by the equivalence of 
(unperturbed) energy and the variable x. 

The resonance state provided by the Lax-Phillips theory in the "energy" representation 
actually therefore corresponds to a distribution of x values in the resonant state. 



° u (x — Re zo) 2 + \Im zo\ 2 ' 

a Cauchy distribution with width \Imzo\. The Cauchy distribution is centered on Rezo, 
corresponding to a shift away from the mean value of x = in the bound state in the 
absence of electric field; as we have seen, the pole moves further to the left with increasing 
field, so that the center of the wave packet moves to the left. 
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